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PHYSICAL ENVIRONMENT OF MASSIVE STAR-FORMING REGION W42 

L. K. Dewangan\ A. Luna\ D. K. Ojha^, B. G. Anandarao®, K. K. Mallick^, 

AND Y. D. Mayya^ 

ABSTRACT 

We present an analysis of multi-wavelength observations from various datasets and Galactic plane 
surveys to study the star formation process in the W42 complex. A bipolar appearance of W42 
complex is evident due to the ionizing feedback from the 05-06 type star in a medium that is highly 
inhomogeneous. The VLT/NAGO adaptive-optics K and L' images (resolutions ~0'.'2-0"l) resolved 
this ionizing source into multiple point-like sources below ^5000 AU scale. The position angle ~15° 
of W42 molecular cloud is consistent with the H-band starlight mean polarization angle which in turn 
is close to the Galactic magnetic field, suggesting the influence of Galactic field on the evolution of the 
W42 molecular cloud. Herschel sub-millimeter data analysis reveals three clumps located along the 
waist axis of the bipolar nebula, with the peak column densities of ^3-5 x 10^^ cm“^ corresponding to 
visual extinctions of Ay ~32-53.5 mag. The Herschel temperature map traces a temperature gradient 
in W42, revealing regions of 20 K, 25 K, and 30-36 K. Herschel maps reveal embedded filaments 
(length ^1-3 pc) which appear to be radially pointed to the denser clump associated with the 05-06 
star, forming a hub-filament system. 512 candidate young stellar objects (YSOs) are identified in 
the complex, ~40% of which are present in clusters distributed mainly within the molecular cloud 
including the Herschel filaments. Our datasets suggest that the YSO clusters including the massive 
stars are located at the junction of the filaments, similar to those seen in Rosette Molecular Gloud. 

Subject headings: dust, extinction - H ii regions - ISM: clouds - ISM: individual objects (W42) - stars: 
formation - stars: pre-main sequence 


1. INTRODUCTION 

Active star-forming regions in molecular clouds often 
contain young star clusters, filaments, bubble(s), and 
massive star(s). Such regions are very promising sites 
for investigation to understand the formation and evolu¬ 
tion of stellar clusters and the interaction between parent 
molecular cloud and embedded massive stars. 

W42 is known as an obscured Galactic gian t H ii region 


et al. 1985) and contains the IRAS 18355—065U source. 


(G25.38—0 .18) towards the inner Galaxy (e.g. Woodward 


Blum et a 


(2000) reported a foreground extinction of 




_ 10 mag m the direction of W42. [Woodward et al.| 

(1985) found a bipolar H ii region with an angular size of 
~12'.'8 (~0.24 pc at a distance of 3.8 kpc) and classified 
its ionizing source as an 07 star. The G25.38—0.18 Hii 
region was further characterized as a core-halo structur e 
with an angular size of ^3' (^3.3 pc) ([Garay et al.|1993|). 
Recently, Spitzer ima.ges (3.6-24 ^m) revealed that W42 
has a bipolar app earance on a much lar ger scale (see 
Figure 10 given in Deharveng et al. 2010), e xtending to 
3(87 (^4.3 pc), and is also known JN39 (e.g., Beaumont 


fc Williams 2010[[Ghurchwell et al.|2006[|Ueharveng et al. 

2010|). The entire complex, which includes an extended 

bipolar nebula associated with G25.38—0.18 H ii region, 
is r eferre d to as W42 complex in this work. [Blum et] 
^ (|2000|) investigated a small embedded stellar cluster 
at the heart of W42 using high spatial resolution near- 
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infrared (NIR) images (hereafter, NIR cluster). They 
studied the NIR K-band spectra of three of the brightest 
stars in the cluster and one of the sources was identified 
as an 05-06 star. It is thought that the W42 complex 
is most likely ionized by this star alone. However, the 
inner circumstellar environment of this source has not 
been studied. The ^^CO profile along the line of sight 
to W42 shows two well-separated velocity components; 
one at 58-69 km s“^ that is physically associated with 
W42, and the other at 8 8-109 km s~^, that is associ¬ 
ated with G25.4NW (e.g. Ai et al. 2013) located at the 


northwest edge of th e extended bipolar nebula (e.g., De¬ 


harveng et al. 2010). Different velocity values suggest 
that these regions are not physically associated a nd are 


no t part of the same star-forming complex (e.g., Ai et 
al. 2013). Using the Cii and ^He radio recombination 


mies, [Cjoheza et al. (2006) measured the velocity of the 
ionized gas to be ~59.6 km s“^ in W42. In the W42 
complex, the velocities of ^^GO gas and the ionized gas 
are consistent with a vel ocity of ^^61 km s~ ^ obtained 


by ammo nia gas (NH3) (jWienen et al. 


al. (2000 


2012 ). 


Blum et| 


reported a distance of ^2.2 kpc to W42, as- 
rat the 0-star is on th e zero-age main sequence 
Anderson et al. ([2009 ) studied the properties of 


summg t. 

(ZAMS). 

molecular cloud associated with the Galactic H ii regions 
including W42 (referred as U25.38—0.18 in their catalog). 
They assigned the a line center velocity of ~64.61 
km s“^ that corresponds to a distance of ~3.8 kpc. The 
near kinematic dis tance reported in the literature (e.g. 


Lester et al. 1985) agrees well with this value. In this 
work, we adopt a distance of 3.8 kpc to the entire W42 
complex. However, we show that our final conclusions 
are not affected by the distance discrepancy. 

Using 870 ^m and 8 fim images of W42 complex, |De-| 
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harveng et al. (2010) reported filaments or sheet-like fea- 
tures along the waist of the bipolar nebula. Additionally, 
they found two dust condensations in the 870 fj,m image 
(in the northern and southern parts) along the waist of 
the bipolar nebula, and one of them (the southern part) 
was se en as an infrared dark c loud (IRDC) in the 8 /im 
image (Deharveng et al.||2010). Methanol maser was de¬ 
tected at 6.7-GJiz in the iiorthern condensatiori . with a 
velocity of ~58.1 km s“^ ( Szymczak et al.||2012). Based 
on the expa nsion of H ii region in a filament, |Deharveng| 


et al. (2010) highlighted physical scenarios to explain the 
existence or the bipolar morphology in N39. 

These previous studies on W42, in general, reveal the 
presence of ongoing star formation, the massive 05-06 
type star and of filament-like features. The study of fila¬ 
ments, their role in star formation process, and the inter¬ 
action of a massive star with its natal molecular cloud are 
yet to be explored observationally in the W42 complex. 
A knowledge of the physical environment of W42 at small 
and large scales is very essential to probe the ongoing 
physical mechanisms. To understand the physical condi¬ 
tions in the complex, we have utilized multi-wavelength 
data covering radio through NIR wavelengths from nu¬ 
merous surveys (e.g. the Multi-Array Galac tic Plane 
Imaging Survey (MAGPIS; Helfand et al.|2006 ), the Co¬ 
ordinated Radio and In frared Survey tor H igliTMass Star 
Formation (CORN ISH; Hoare et aL| 2012|), the Galactic 
Ring Survey (GRS;|jackson et al.|20(J6p, the APEX Tele- 
scope Large Area Survey ot the Galaxy (ATLASGAL; 
Schuller et al.|200 9[), the fferscfeel Infrar ed Galactic Plane 
i-GAL, Molinari et al. 20R) ), the MIPS Inne r 


Survey (Hi- 

Galactic Plane Survey (MIPSGAL; |0arey et al.||2005|) 
the Galactic Legacy Infrared Mid-Plane Su rvey Extraor- 
dinaire (GLIMPSE; [Benjamin et al. 2003), the UKIRT 


et al. 


Wide-held I nfrared Survey tor H2 (UWISH2; Froebrich 
201 1|), the UKIRT NIR Galactic Plane Survey 
Lawrence et al.|2007|) , the Galactic Plane Infrared 




Polarization Survey (GpIPS; Clemens et al. 2012), the 
ESO Very Large Telescope (VL'L) archi ve, and the Two 


Micro n All Sky Survey data (2MASS; Skrutskie et al. 


2006)). The ESO-VLT archival adaptivAoptics INIR im- 


ages are used to trace the small scale environment (inner 
5000 AU) of the 05-06 star. All these surveys are uti¬ 
lized to explore the different components associated with 
the complex, viz., molecular gas, ionized emission, dust 
(warm and cold) emission, shocked emission, distribution 
of dust temperature, column density, extinction, mag¬ 
netic held, and embedded young stellar objects (YSOs). 

This paper is organized as follows. In Section we 
describe in detail various datasets along with reduction 
procedures. Section focuses on the results related to 
the physical environment and point-like sources (classi- 
hcation and their analysis) of the region from various 
datasets. The possible star formation scenarios are dis¬ 
cussed in Section Gonclusions are presented in Sec¬ 
tion [5] 


2. DATA AND ANALYSIS 

We have used multi-wavelength data to study the star 
formation and feedback of massive star on the surround¬ 
ing interstellar medium (ISM) in the W42 complex. The 
size of the selected area is ~ 1514 x 1414, centered at 02000 
= 18^38'"14®, ^2000 = —06°47'24", which corresponds to 
a physical scale of ~17 pc x 16 pc. 


2.1. NIR Data 

NIR photometric data were obtained from the UKIDSS 
6‘^ archival data release (UKIDSSDR6plus) of the GPS. 
UKIDSS observations were per formed using the U KIRT 
Wide Field Gamera (WFGAM; |Casali et al.||2007|). The 
hnal huxes were calibrated using the 2MASS data. A full 
description of data reduction and calibration pro c edure s 
are given in Dye et al. (2006) and Hodgkin et al. (2009), 
respectively. In order to select only reliable point sources 
fro m the GPS c a talog, we utilized the criteria suggested 
by Lucas et al. (2008). For the sources common in all 
the three JMR (jR A) bands and those detected only in 
the H and A bands, separate criterigQ were adopted for 
selection. These criteria include the removal of saturated 
sources, non-stellar sources, and unreliable sources near 
the sensitivity limits. Following the criteria, our selected 
GPS catalog contains sources fainter than J = 12.5, H = 
11.6, and K = 10.2 mag to avoid saturation. The magni¬ 
tudes of saturated bright sources were retrieved from the 
2MASS catalog. We selected only those sources which 
have magnitude error of 0.1 or less in each band, to ob¬ 
tain good photometric quality. Following the above cri¬ 
teria, we found 17745 sources common to all the three J, 
H, and K bands. In addition to these detections, 8754 
sources having only the H and K bands detection were 
also obtained. 

We retrieved archival adaptive-optics imaging data 
towards W42 region from the ESO-Science Archive 
Facility (ESO proposal ID: 089.C-0455(A); PI: Joao 
Alves). Observations were made with 8.2m VLT 
with NAOS-GONIC A (NAGO) adap tive-optics system 
(Lenzen et al. 20031 Rousset et al. 2003) in Kg-band 
(Ac = 2.18/rm, AA = 0.35 /im) and L'-band (Ac = 
3.80p,m, AA = 0.62/rm). We obtained five Kg frames 
and six L' frames of 24 and 21 seconds of exposures, 
respectively. The final processed NAGO images were ob¬ 
tained through the standard analysis procedure such as 
sky subtraction, image registration, combining with me¬ 
dian method, and astrometric calibration, using IRAF 
and STAR-LINK softwares. The astrometry calibration 
of NAGO images was performed using the GPS K-band 
point sources. NAGO Kg and L' images have plate scales 
of O'.'054/pixel and 07027/pixel, respectively, with reso¬ 
lutions varying from 0'.'2 (-760 AU) - 0"1 (-380 AU). 


2.2. H 2 Narrow-band Image 

Narrow-band H2 (v = 1 — 0 S(l)) image at 2.12 fim 
was retrieved from the UWISH2 database. To obtain a 
continuum-subtracted H2 map, the point spread function 
of GPS K-band image was matched and scaled to the H2 
image. 


2.3. H-hand Polarimetry 

Imaging polarimetry in NIR H-band (1.6 /rm) was ob¬ 
tained from the GPIP^ Observations were taken with 
the Mimir instrument, on the 1.8 m Perkins t elescope, 


in H-band l inear imaging polarimetry mode (see Glemens 
et al.|[20T2 for more details). The polarization data to¬ 


wards W42 complex were covered in two GPIPS fields 


^ In the appendix jA| we give the SQL script that is especially 
written to implement these criteria. 

^ http://gpips0.bu.edu/ 
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i.e., GP0612 {I = 25°.319, b = -0°.240) and GP0626 {I 
= 25°.447, b = —0°.165). Only those sources with Usage 
Flag (UF) = 1 and P/cTp > 2 were selected for the study, 
to obtain good polarimetric quality, where P is the polar¬ 
ization percentage and ap is the polarimetric uncertainty. 
These conditions provided a total of 234 stars from two 
GPIPS fields. 

2.4. Spitzer Data 


Spi tzer Infrared Array Camera (IRAC; Fazio et al. 
20041 3.6-8.0 /im images were obtained from the 
GLIMPSE survey. We obtained the photometry of point 
sources from the GLIMPSE-I Spring ’07 highly reliable 
catalog. The GLIMPSE-I catalog does not provide pho¬ 
tometric magnitudes of some sources, which are well de¬ 
tected in the images. Aperture photometry was per¬ 
formed for such sources using the GLIMPSE images at a 
plate scale of Cf.'6/pixel. The photometry was extracted 
using a 274 aperture radius and a sky annulus from 274 to 
773 in IRAPq Apparent magnitudes were calibrated us¬ 
ing the IRAG zero-magnitudes (i.e., 18.59 (3.6 /rm), 18.09 
(4.5 /.tm), 17.49 (5.8 fim), and 16.70 (8.0 pTo)) including 
aperture corrections (see IRAC Instr ument Handbook 


(Vers ion 1.0, 2010 Eebruary) and also Dewangan et al. 
20121. We also utilized MIPSGAL 24 ^m photometry 


in this work. We performed aperture photometry on the 
MIPSGAL 24 pm image to extract point sources. The 
photometry was obtained using a 1" apert ure radius and 


a sky ann ulus from 7" to 13" in IRAF (e.g. Dewangan et, 
ar]|2015a). MIPS zero-magnitude flux density, including 


aperture correction was used for the photometric cali¬ 
bration, as reported in the MIPS Instrument Handbook- 
Ver-fl 

2.5. Herschel and ATLASGAL Data 

Herschel Hi-GAL continuum maps were retrieved for 
70 pm, 160 pm, 250 pm, 350 pm, and 500 pm. The 
beam sizes of these bands are 578, 12", 18" , 25", and 37" 
(Poglitsch et al. 2010 Griffin et al. 2010), respectively. 
We obtained the processed level2_5 products, using the 


Hersc hel Interactive Processing Environment (HIPE, Ott 
2010). The images at 70-160 pm were in the units of Jy 
pixel” ^, while the images at 250-500 pm were calibrated 
in the surface brightness unit of MJy sr”^. The plate 
scales of the 70, 160, 250, 350, and 500 pm images are 
3.2, 6.4, 6, 10, and 14 arcsec/pixel, respectively. 

The sub-millimeter (mm) continuum map at 870 pm 
(beam size ~1972) was obtained from the APEX AT¬ 
LASGAL archival survey. 

2 .6. (J=l—0) Line Data 

The (J=l—0) line data were downloaded from 

the GR^ The survey data have a velocity resolution 
of 0.21 kms”^, an angular resolution of 45" with 22" 
sampling, a typical rms sensitivity (Icr) of « 0.13 K, a 
velocity coverage of —5 to 135 km s~^, and a m ain beam 
efficiency (r/mb) of ~0.48 ( Jackson et al.|[2006 |. 


2.7. Radio Centimeter Continuum Map 

Radio continuum map at 20 cm (beam size ^6") was 
obtained from the VLA MAGPl^ CORNlSt0 5 GHz 
(6 cm) high-resolution radio continuum data (beam size 
^175) are also utilized in this work. CORNISH 5 GHz 
compact source catalog was also retrieved from |Purcell| 
eFaLlllM^. 


2 .8. Other Data 

We used the observed positions o f the 6.7-GHz 
meth anol maser (iSzymc zak et al. 2012) and the 05-06 
star dBlum et al.||2000|). 'Lhe ^operties of th e molec¬ 
ular cloud U25.38—0.18 (Anderson et al. 2009) are also 
utilized in this work. 

3. RESULTS 

3.1. Multi-phase environment of W42 complex 
3.1.1. The bipolar nebula and its heart 

The spatial distribution of warm dust emission traced 
in Herschel 70 pm image towards our selected region 
(size ~ 15(4 X 1414) is shown in Figure]^ which depicts 
the previously known bipolar nebular morphology. The 
bubble region appears to extend several parsecs (~11 pc 
X 7 pc). As traced in 8.0 pm emission, the waist and the 
edges of the nebula are highlighted in Figure (also see 
Figure 10 in Deharveng et al. 20T0|. The 24 pm image 
also traces the warm dust emission and is saturated near 
the location of the 05 -06 star (see Figure 10 given in 


Deharveng et al. 2010). Furthermore, the 24 pm emis- 
sion is enclosed by the 8.0 pm emission. The G25.4NW 
region lies in the direction of the W42 complex, but is not 
associate d with it, as ca n be seen from its very different 
velocity (Ai et al.||2013 ) which indicates a very different 
distance. 'Lheretore, the results of G25.4NW region are 
not discussed in this work. The radio continuum emission 
at MAGPIS 20 cm (beam ^^6"), which traces the ionized 
emission, shows the central ~4.28 pc part of the bipolar 
nebula. The peak of 20 cm emission is found to be ap¬ 
proximately coincident with a well-classified 05-06 star 
(see Figure [^. This source could be the powering source 
of the entire extende d emission. Using the MAGP IS 20 
cm continuum data, Beaumont & Williams (2010) esti¬ 
mated the Lyman continuum photon number (iogN,^^) 
within the N39 nebula to be ^49.24 for a distance of 3.7 
kpc (or logN„„ ^49.26 for a distance of 3.8 kpc). This 
observed value is in agreement with the theoretical 
value for a spect ral type of 05 star (logN„„ ^49.26; 


Martins et al. 2005). In this case, we assume that no 
ionizing photons are absorbed by the dust in the ionized 
gas, which is probably not true because of the presence of 
warm d ust emission i nside the Hii region. As mentioned 
before, Blum et al. (2000) have reported a distance of 
2.2 kpc to W42. We find logN,j„ ~48.79 for a distance 
of 2.2 kpc, which corresponds to a single ionizing star of 
spectral type 06.5V (see Table 1 in|Martins et al.|(|2005) 
for theoretical values). The radio spectral type of the 
exciting star at 2.2 and 3.8 kp c is consistent wit h a MK 
type spectrum of 05-06 (e.g. Blum et^^|2000 ). 

The inset on the top right of I'igure shows the heart 
^ IRAF is distributed by the National Optical Astronomy Ob- of the bipolar nebula, aS Seen in the opiteer-IRAG 5.8 
servatory, USA 

^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/mips/mipsinstrumei^tBidt^l^i^hki/f’d.ucllnl.org/gps/index.html 
^ http://www.bu.edu/galacticring/ ^ http://cornish.leeds.ac.uk/public/index.php 
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/im image overlaid with MAGPIS 20 cm and IRAC 5.8 
fim contours. The 5.8 /rm contour emission reveals a 
structure (a spatial extension of ~1 pc scale), w hich was 
interpreted as an i onized cavity-like structure by |Dewan-| 
gan et al.| (|2015b|). This structure is well traced in H 2 
2.12 fj,m emission and i n the continuum emission at 3.6- 
8.0 fim, and 20 cm (see Dewangan et al.||2015b|. COR¬ 
NISH 5 GHz high-resoliition radio cohtmuum emission 
(beam ~1"5; Purcell et al. |2013|) traces two compact 
radio sources (i.e. GU25.3824—00.1812 (angular scale 

^3"1) and G025.3809—00.1815 (angular scale ~8"3)) lo¬ 
cated inside this cavity-like structure (see Figure |^). 
CORNISH 5 GHz radio contours are overlaid on the 
GPS K band image in Figure!^. A comparison of ra¬ 
dio emission between MAGPIS20 cm (beam ^6") and 
CORNIS H 5 GHz (beam ~I"5) can be seen in Figures]^ 
and [ 2 ^. Purcell et al. (20131 also measured the inte¬ 
grate flux densities equal to 200.13 mJy and 460.83 mJy 
for G025.3824-00.1812 and G025.3809-00.1815 sources, 
respectively (see Table for coordinates). The spec¬ 
troscopically identified ^5-06 star is located within 
the extension of G025.3809—00.1815. The peak posi¬ 
tions of these radio sources, G025.3824—00.1812 and 
G025.3809—00.1815, are ^5"4 and ~5"5 away from the 
location of the 05-06 star, respectively. In order to infer 
the spectral class of each of the compact radio sources, 
we have estimated the number of Lyman continuum pho¬ 
tons (N„„) using the integrated flu x density following the 


equation of Matsakis et al. (19761 (also see Dewangan et 
al. 2015a for more details). I'he calculations were car- 


ried out for a distance of 3.8 kpc and for the electron 
temperature of 10000 K. We compute N„„ (or logN„„) to 
be -2.6 X 1047 s-4 (47.41) and -5.9 x lO^^ s"! (47.77) 
for G025.3824-00.1812 and G025.3809-00.1815, respec- 


tively. These values correspond to a single ionizing star of 

spectral type BOV (see Table 1 in Smith et al 

1 (12002) for 

theoretical values) and 09.5 V (see 'Fable 1 in 

Martins et 

al. 

anc 

(20051 for theoretical values) for G025.382 
GU25.3809—00.1815, respectively. These i 

4-00.1812 

adio peaks 


are individual Hii regions. Note that the extended ra¬ 
dio emission at 20 cm does not appear to be produced 
by these BOV and 09.5V type stars together. There 
is no infrared counterpart to the peak of radio source 
G025.3809—00.1815 seen in the GPS NIR image, while a 
NIR counterpart of G025.3824—00.1812 is found in the 
GPS NIR image. Note that in the presence of the 05- 
06 star, there will hardly be any noticeable additional 
effect of BOV spectral type star, because the radio flux 
of BOV type star (—0.2 Jy) is about 1/80*^ of the 05-06 
type star (i.e. —16 Jy). Therefore, the W42 complex ap¬ 
pears more likely to be excited by the spectroscopically 
identified 05-06 type star. 

In order to study the small scale environment of the 
05-06 type star, we utilized the ESO-VLT archival 
NIR images. In Figure [^, we present the VLT/NAGO 
adaptive-optics NIR images towards G025.3809—00.1815 
(including the 05-06 star) in Kg and L' bands. Fig¬ 
ure shows the VLT/NAGO Kg image towards 
0025:3824-00.1812 radio source. The VLT/NAGO im¬ 
ages provide the circumstellar view of the 05-06 star 
in inner regions within 5000 AU. Within this scale, the 
05-06 star is resolved into at least two stellar sources 
in the Kg image, while the L' image resolves the 05-06 


star into at least three stellar sources. We do not find 
any nebular feature in NAGO images within 10000 AU of 
the 05-06 star. In general, it is observationally known 
that mas sive stars are often found in binary and multiple 
systems (Duchene & Kraus 2013|) . In M8 massive star¬ 
forming region, jOoto et al.| (|2U06|) resolved an “O” star 
(designated as ^"Herschel — 36”) into multiple sources 
using the VLT/NAGO images. It appears that the 05- 
06 star in W42 is associated with a very similar system 
of multiple stellar sources lik e Herschel — 36 in M8 (see 


Figure I in Goto et al. 2006). The VLT/NAGO Kg im¬ 
age shows at least 5 stellar sources towards the peak po¬ 
sition of G025.3824—00.1812, within a scale of —10000 
AU. There are no nebular features seen towards these 
sources in NAGO images. We do not have detections 
of these resolved sources in both NAGO bands, hence 
we cannot provide any quantitative predictions (s uch as 
color a nd sp ectral type) in this work. Previously, Blum 
et al. (2000) detected only three point-like sources to¬ 
wards G025.3824— 00.1812 and studied the K-band spec¬ 
tra of one of the three source s (see source #3 marked in 
Figure 1 in Blum et al.|2000|). These authors did not de¬ 


tect any stellar absorption features, and also found NIR 
excess emission associated with source, which led them 
to suggest this source as a YSO candidate having ex¬ 
tinction of about 10 mag. High resolution spectroscopic 
study will be helpful to identify main ionizing source of 
the complex. On the other hand, it is also probable that 
the W42 complex is being ionized by a small cluster of 
O and early B type stars rather than by a single star. 
Additionally, Dewangan et al. (2015b) reported a parsec 
scale H 2 outflow that is driven by an infrared counterpart 
of the 6.7-GHz methanol maser emission, namely, W42- 
MME. W42-MME is observed at wavelengths longer than 
2.2 ^m and is classified as a deeply embedded massive 
YSO (stellar mass —1 9±4 and extin ction —48±15 
mag) (see Figure (|Dewangan et al.| |2015b|) . They 
also investigated a jet-like feature in W42-MME using 
the VLT NIR adaptive-optics images. The 05-06 star is 
located at a projected linear separation of about 0.22 pc 
from the W42-MME. These observational features sug¬ 
gest that massive star formation is cur rent ly taking place 
in the W42 complex (also see Section 4.3). 


3.1.2. Tracers: warm dust, cold dust, and molecular gas 

In Figure we show a longer wavelength view of the 
W42 complex. The images are shown at 160 /rm, 250 /im, 
350 /rm, 500 /rm, 870 //m, and integrated ^^CO (J=l—0) 
intensity contour map. The CO map is integrated in the 
[58,69] km s“4 velocity range. The 160- 870 /rm emission 
traces cold dust components (see Section 3.2 for quantita¬ 
tive information). Most of the emission m all the maps 
is concentrated toward the waist of the bipolar nebula 
(i.e., two condensations, namely northern and sou thern), 
as previously reported by Deharveng et al. (2010). Both 
of these condensations have a similar position angle of 
—15°. CO data are very useful to study the morphology 
of molecular cloud. In Figure [Sf, an elongated molecular 
cloud structure is evident along the waist of the neb¬ 
ula, where the northern and southern condensations are 
clearly visible as prominent bright features. The north¬ 
ern dust component appears to contain two CO clumps 
(cNl and cN2) and is associated with previously known 
NIR cluster including the 6.7-GHz maser and the 05-06 
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type star (see Table for coordinates). The southern 
condensation, traced as a CO clump (cSl), is found to 
be linked with the IRDC seen at 8.0 ^m. The IRDC ap¬ 
pears as a bright filament (length ~5 pc) in emission at 
wavelengths longer than 70 ^m. In addition to the north¬ 
ern and southern condensations, the filamentary features 
(labels: a, b, c, d, e, f, and h) are also present in the 
Hi-GAL maps. In the proximity of the northern conden¬ 
sation, the filamentary features are highlighted based on 
visual inspection of the Hi-GAL 250 /im map (see curves 
in Figure]^; hereafter Herschel filaments). The posi¬ 
tion angles of these filaments (~40-170 degrees) are also 
shown in Figure]^. The lengths of these filaments vary 
between 1 and 3 pc scales. The orientations of these 
filaments are different from the position angles of the 
northern and southern condensations. Note that these 
Herschel filamentary features are not detected in the 870 
/rm map (see Figure [^). 

Taken together. Figures [T}|^ illustrate the location of 
the 05-06 star, 5 GHz radio peaks, 6.7-GHz maser, ion¬ 
ized emission, molecular emission, filaments, IRDC, and 
dust (warm and cold) emission in the complex. 


3.1.3. H 2 emission 

In Figure]^, we display a continuum-subtracted 2.12 
/rm H 2 (v=I — 0S(I)) image near the waist of the neb¬ 
ula (size of the selected region ~5.7 pc x 5.4 pc). The 
map shows H 2 emission toward the waist and edges of 
the bipolar nebula. The comparison of H 2 emission with 
the 8.0 /rm features can be seen in Figures and 
The 8.0 /xm emission, which contains 7.7 /xm and 8.6 
/xm polycyclic aromatic hydrocarbon (PAH) features in¬ 
cluding the continuum, traces a photodissociation region 
(PDR). Considering the distribution of H 2 emission, the 
ionized gas (see Figure 0 and the 8.0 ^m emission, we 
suggest that the H 2 emission likely traces the PDR in 
the complex. However, it is also possible that the H 2 
emission probably originates in t he s hock due to the ex¬ 
panding Hii region (see Section 3.3). The origin of the 
H 2 emission in terms of shocks and UV fluorescence is 
often explained in the literature by the observed ratio 
of the l-OS(I) to the 2-IS(I) intensity. However, we do 
not have observations of the source in the 2-IS(I) line 
of H 2 filter. Therefore, we cannot conclude the origin of 
H 2 emission in this work. Additionally, the H 2 emission 
near the 6.7-GHz maser emission at a parsec scale can¬ 
not be ruled out due to outflow activi ty (see Dewangan 
et ah 20I5b[ and references therein). Lee et ai.l (2U14|) 
and|Sninn et al.| (|20I4|) reported the detection of [h'e H| 
emission in the southwest of the 6.7-GHz maser emission, 
which fur ther suggests the presence of shocks (see Fig¬ 


ure 12 in Lee et al. 


[ 2 ^. Note that the H 2 emission 


is also detected at the tip of the southern condensation 
(see arrow in Figure]^). The 8.0 /xm and the H 2 emis¬ 
sions have very similar spatial structure at the tip of the 
southern condensation. 

In summary, the presence of H 2 features in the com¬ 
plex provides the observational signatures of the outflow 
activity as well as the impact of the UV photons. 


3.2. Temperature and column density maps of WfE 
complex 

In previous sections, we qualitatively studied the mor¬ 
phology of the complex using Herschel images. In this 


section we present the temperature and column density 
maps of the complex, derived using Herschel images. A 
knowledge of the temperature distribution is important 
to infer the physical conditions in the cloud. The distri¬ 
bution of column density allows to infer the extinction, 
mass, and density variations in the cloud. The final tem¬ 
perature and column density maps were o btained follow¬ 
ing th e same procedures as described in [Mallick et al. 


(2015). These maps were determined from a pixel-by¬ 
pixel spectral energy distribution (SED) fit with a mod¬ 
ified blackbody curve to the cold dust emission in the 
Herschel 160-500 /xm wavelengths regime. We did not 
include Herschel 70 /xm data, because the 70 /xm emis¬ 
sion comes from UV-heated warm dust. 

Here, we give a brief description of the procedures. In 
the first step, we converted the surface brightness unit 
of 250-500 /xm images to Jy pixeU^, same as the unit of 
160 /xm image. Next, the final processed 160, 250, and 
350 /xm images were convolved to the angular resolution 
of the 500 /xm image (~37"), using the convolution ker¬ 
nels available in HIPE software and then regridded on a 
14" raster. We then estimated a background flux level. 
The dark and featureless area far from the main cloud 
complex is selected for the background estimation. The 
background flux level was obtained to be -2.98, 1.14, 0.55, 
and 0.19 Jy pixel“^ for the 160, 250, 350, and 500 /xm im¬ 
ages (size of the selected region ~ 6.'6 x 713; central coor¬ 
dinates: aj 2 ooo = 18H2'"34L6, 6 j 2 ooo = -06°23'11".4), 
respectively. Finally, a modified blackbody was fitted to 
the observed fluxes on a pixel-by-pixel ba sis to generate 


the m aps (see equations 8 and 9 given in Mallick et al. 


2015). The fitting was done using the four data points 


tor each pixel, keeping the dust temperature (T^^) and 
the column density (V(H 2 )) as free parameters. In the 
calculations, we adopted the me an molecular weight per 
hydrogen molecule {fJ,H 2 =) 2.8 ( Kauffmann et al.]2008[ 
and an absorption coefficient {ki, =) 0.1 (xx/1000 GHz) 
cm^ g“^, including a gas-to-dust ra tio (Rf =) of 100, with 
a dust spectral index of ,0 = 2 (see Hildebrand||1983 ). 

The final temperature and column density maps (m- 
gular resolution ~37") are shown in Figures andj^, 
respectively. The temperature map clearly Snows tem¬ 
perature variations in the complex. Our temperature 
map shows considerably warmer gas (T^; ~30-36 K) to¬ 
wards the heart of the bipolar nebula, where the ionizing 
source is located. We find that the mean tem perature of 
mo st of t he nebula is ~25 K in the complex. |Wienen "et 


al. ( 2012 ) derived the line parameters (such as the gas ki- 


netic temperature and rotational temperature) towards 
the clumps found by the ATLASGAL survey. We found 
one ATLASGAL dust clump in our selected re gion (see 
Figur e j^), w hich has NH 3 line parameters from Wienen[ 


et al. (2012). The rotational temperature of this clump 


was found to be ^20 K. It is consistent with our esti¬ 
mated temperature value towards this clump. G25.4NW 
region also dominates in the temperature map with a 
peak temperature of ^45 K. Two clumpy as marked in 
the CO map (cNl and cN2; see Figure [^), are traced 
in the northern condensation with peak column densi¬ 
ties of ~3.6 X 10^^ and 3.0 x 10^^ cm“^, which cor¬ 
respond to visual extinctions of A y ^38.5-32 nia.g, as- 
suming the classical relation from Bohlin et al. ( 2000| 
(i.e. Ay = 1.07 X 10“^^ V(H 2 )). The colder gas (T^ 
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^20 K) is found near the column density peak (5 x 10^^ 
cm“^; Ay ^53.5 mag) in the southern condensation 
clump (cS l), where the IRDC is s een in 8.0 /xm map (see 
Figure f^. | Anderson et al.| (|2009|) studied the molecular 
cloud 1125.38—0.18 (i.e. northern condensation) using 
the GRS CO data and estimated its column density to 
be ~1.05 X 10^^ cm“^, which is well in agreement with 
our column density value towards the northern condensa¬ 
tion. The condensation associated with G25.4NW region 
has the highest column density, as traced in Figure [^. 
Herschel filaments are also drawn in Figure Due to 
the coarse resolution of the column density map, some of 
the filaments are not resolved as seen in the 250 /xm map. 
The value of column densities towards these filaments is 
found to be ~1.5 x 10^^ cm“^. The column density map 
shows the dense clump (cNl; peak fV(H 2 ) ~3.0 x 10^^ 
cm“^) associated with the 05-06 star and W42-MME, 
where several filaments (“a-h”; see Figure]^) appear to 
be radially directed to this clump (see Figure]^ and also 
Figure!^, re vealing a “hub-filament” morphology (e.g. 


Myers P09D . We should mention here that all these fil¬ 


aments areTocated well within the W42 molecular cloud, 
except the two filaments “h” and “g” that appear to be 
extended towards G25.4NW from the northern conden¬ 
sations (also see Figurej^. We suggest that this could be 
a projection effect because these two regions (W42 and 
G25.4NW) are not physically linked. Four filaments (e, 
f, g, and h) are also seen in the H 2 and 8.0 /xm maps, 
while remaining filaments (a, b, c, and d) are traced 
only in the Herschel images. These filaments are nei¬ 
ther perfectly parallel nor perpendicular to the parental 
molecular cloud. The cavity-like structure (embedded 
in cNl clump) appears at the junction of filaments (see 
Figure Jm. This cavity seems to be an important link 
for uncferstanding the ongoing physical processes in the 
complex. 

Using the column density map (in Figure [^), we 
also estimated the total masses of the clumps associ¬ 
ated with the northern condensation (cNl and cN2) and 
the southern condensation (cSl). Total column densi¬ 
ties (EAf(H 2 )) of these clumps w ere computed us i ng th e 


“GLUMPFIND” IDL program (Williams et al. 
Following t he relation M oc Eiv (142) (also see 


1994). 


Vlallick 


et al. 2015), the masses of the clumps cNl, dNI2, and 
cSl are computed to be ~730, ~840, and ^2564 M©, 
respectively. 


3.3. (J=l—0) kinematics in W42 complex 

The profile showed that the W42 complex is cov¬ 
ered in the velocity range of 58-69 km s“^. In Figure 
we present the integrated GRS (J=l—0) velocity 

channel maps (at intervals of 1 km s“^), which reveal the 
morphology and different molecular components along 
the line of sight. The velocity channel maps trace the 
prominent northern and southern condensations in the 
complex. The ionizing source is located within the north¬ 
ern condensation. The physical association of the ioniz¬ 
ing source was confirmed by the velocities of molecular 
and ionized gas, as mentioned in the introduction. Addi¬ 
tionally, the maps reveal the regions empty of molecular 
GO gas (see highlighted regions in Figure]^. 

The integrated ^^GO intensity map and position- 
velocity maps are shown in Figure The position- 


velocity maps (right ascension-velocity and declination- 
velocity) show the presence of a noticeable velocity gra¬ 
dient in both the northern and southern molecular com¬ 
ponents. The position-velocity plots of ^^GO gas re¬ 
veal an almost semi-ring-like or inverted C-like structure 
(see Right Ascension-velocity panel in Figure]^. Such a 
structure in the position-velocity plot is consisten t with 
the model fo r an expanding shell (Arce et al.|2011|) . Arce 
et al.| (|2011|) performed modeling of expanding bubbles 
m a turbulent medium and compared with the observed 
structures in Perseus molecular cloud. They suggested 
that the semi-ring-like or G-like structure is characteris¬ 
tic of an expanding shell. These authors also pointed out 
that a ring-like structure can be evident in the position- 
velocity plot when the powering source is located at the 
center of the region. Since the complex harbors a power¬ 
ful 05-06 star. Therefore, the existence of an inverted 
G-like structure can be explained by the expanding H ii 
region. Following the position-velocity maps, we infer 
the expansion velocity of the gas to be ~3 km s“^. 

We further analyzed the gas distribution in the north¬ 
ern component using the position-velocity analysis and 
found the receding gas (65-68 kms“^), approaching gas 
(59-63 kms“^), and rest gas (63-65 kms“^) components 
(see bottom left panel in Figure]^. Due to the coarse 
beam of GO data (beam size ~45^, we cannot pinpoint 
the exact exciting source of this outflow. This outflow 
signature is associated with whole northern condensa¬ 
tion, hence it could be related to the expanding Hii re¬ 
gion associated with the 05-06 star. 

In order to trace out the regions of direct interaction of 
ionized gas, we carefully analyzed the position-velocity 
maps (see Figure and found two cavities empty of 
molecular gas in the south-west and south-east directions 
with respect to the 05-06 star (see Figures and [I^. 
These cavities are traced in the velocity ranges of 63- 
68 kms“^ and 58-64 kms“^ in the south-west and south¬ 
east directions, respectively. The spatial locations of the 
cavities show the symmetry with respect to the main 
elongated molecular cloud. The existence of the molec¬ 
ular cavities suggests that the molecular gas has been 
eroded by the ionizing gas. 

Note that the GRS ^^GO data do not allow us to ex¬ 
plore any outflow signatures towards W42-MME due to 
the coarse beam (beam size ~45"). Hence, high reso¬ 
lution molecular line observations are required to obtain 
better insight into the molecular outflows in the complex. 
Gombining the inferences from the GO kinematics, the 
influence of ionizing gas on the surroundings is evident. 


3.4. Dynamical age of the Hii region and feedback of a 
massive star 

In this section, we compute the dynamical age of the 
H II region. Knowledge of this physical parameter helps 
in understanding the local star formation process by the 
interaction of the H ii region with the surrounding ISM. 
Here, we performed the calculations based on the radio 
continuum flux at 20 cm and radio spectral type of the 
ionizing source. The dynamical age (tdyn) of a spherically 
expanding H ii region of a radius R is est i mated using 
the model described by Dyson & Williams (|1980) and is 
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given by: 


^dyn 


7 Cs 


7/4 


- 1 


( 1 ) 


where Cg is the isothermal sound velocity in the ionized 
gas (Cg = 10 km s“^) and Rg (= 3N„^/47rnoa_B)^^^ is 
the radius of the initial Stromgren sphere, where the ini¬ 
tial particle number density of the ambient neutral gas is 
“no”, the radiative recombination coefficient is “ag” (= 


2.6 X 10"^^ X (10-* K/Te)0 '^ cm3 g-i. ggg 
and N„ 


Kwan 1997), 

is the Lyman continuum photons per second. 
In this calculation, we assume that the Hii region as¬ 
sociated with the complex was spherical in morphology 
during its initial phase and, with time, it evolved into 
the homogeneous surrounding environments. Here, we 
use N,„, = 1.83 X 10"^® s”^ or logN„^ = 49.26 (see Sec¬ 
tion 


3.1.1) for a radius of the Hii region, R 


2.14 pc 

(Beaumont & Williams 2010), and as = 2.6 x lO”^^ 
cm3 g-i ^ 10000 K. Anderson et al. (2009) esti¬ 

mated the column density ( 1 N(H 2 )) and size (L) of the 
molecular cloud U25.38—0.18 (i.e. northern condensa¬ 
tion) as ^1.05 X 10^^ cm“^ and 1.'5 (5.12 x 10^® cm 
at a distance of 3.8 kpc), respectively. The H 2 number 
density is computed to be ^2051 cm”® using the relation 
Nh 2 {cm~^)/L (cm), for an assumed spherical structure. 
Using the values of N„„, ng, and R in Equation 1, we ob¬ 
tain Rg and dynamical age of the Hii region as ~0.51 pc 
and ^0.32 Myr, respectively. If we also adopt a distance 
of 2.2 kpc to W42 then we obtain a radius of the Hii 
region (R) « 1.24 pc , N„,„ = 6.13 x 10^® s”® or logN„^ 
= 48.79 (see Section 3.1.1), size of the molecular cloud 
(L) ^2.97 X 10^® cm, H 2 number density ~3541 cm”®, 
Rg ~0.35 pc, and dynamical age of the Hii region ^0.22 
Myr. The dynamical age of the H ii region calculated at 
2.2 kpc is estimated to be about 69% of that derived at a 
distance of 3.8 kpc. The estimated dynamical age should 
be considered with some caution, because of the assump¬ 
tions involving spherical geometry and uniform density 
distribution. Also, the observed bipolar morphology of 
the complex could have originated due to the evolution 
of the Hii region in a medium with strong density gra¬ 
dients. 

In the previous sections, we found that the massive 
star (05-06) is located within the northern conden¬ 
sation. In order to study the feedback of this massive 
star on the southern condensatio n, we calculated th e 
following pressure components (e.g. Bressert et al.|20I2): 
(i) pressure of an H ii region {FinT) = 

w 


liruHC^s 


3Nuv 

477 CKB 


n-, 


(ii) radiation pressure (Prad) = Lboi/^T^cD 
(hi) stellar wind ram pressure (Pwind) = 

(iv) pressure exerted by the self-gravity of the surround¬ 
ing molecular g as {Psdoud) ~ ttGE^ (e.g. Harper-Clark 
&: Murray||2009| ). 

In the relations above, N„j, and as are defined as in 
Equation 1, = 2.37 (approximately 70% H and 28% 

He by mass), m/f is the hydrogen atom mass, Cg is the 
sound speed of the photo-ionized gas (= 10 km s”®), 
Mw is the mass-loss rate, Vu, is the wind velocity of 
the ionizing source, Lhoi is the bolometric luminosity 


of the region, E (= Mgcioud/T^Rl) is the mean mass 
surface density of the southern condensation, Msdoud 
is the mass of the molecular gas associated with the 
southern condensation, and Rg is the radius of the 
molecular region. The pressure components associated 
with massive star are evaluated at Dg = 2.14 pc (radius 
of the Hii region) on the southern condensation from 
the position of the 05-06 type star. The luminosity of 
the exciting 05 type star is ^3.2 x 10® Lfr^ (l ogL/Lfr^ = 
5.51; see Table 1 given in Martins et al. 2005). 

S ubsti tuting Mgdoud (clump cSl) « 2564 JvIq (see Sec¬ 
tion 3.2), Rc « 1.6 pc, Lboi = 3. 2 x 10® L( 7 ), = 2 .0 x 

10” Mq yr”® (for an 06V star; de Jager et al.|l988 ), 

= 2500 km s”® (for an 06V star; Prinja et al. 1990) in 


the above equations, we find a surface density E « 0.067 
g cm”2, Psdoud « 9.3 X 10”®° dynes cm”^, Phii « 9.3 
X 10”®° dynes cm”^, P^ad ~ 7.57 x 10”®® dynes cm”^, 
and Pwind ~ 5.76 x 10”®^ dynes cm”^. The compari¬ 
son of different pressure components associated with the 
massive star suggests that the pressure of the H ii region 
is relatively higher than the radiation pressure and the 
stellar wind pressure. 

Note that a typical cool molecular cloud (temperature 
^20 K and particle density ~10®-10^ cm”®) has pres¬ 
sure values ~10”®®-10”®^_dynes cm”^ (see Table 7.3 


of Dyson & Williams 1980). We estimate the value of 
Psdoud ~ 9.3 X 10”®° dynes cm”^, which is relatively 
higher than the pressure associated with a typical cool 
molecular cloud. It suggests that the surrounding molec¬ 
ular cloud has been compressed to increase the pressure. 

Eurthermore, we also derived the virial mass of the 
southern condensation using the line width (Av) of ®®CO 
velocity profile obtained usi ng a single Gaussian fitt ing. 
The virial mass is given in MacLaren et al. ( |1988 ) as 
Myij. (M 0 )=kRcAv^, with Rc (in pc) is the radius 
of the clump as defined above, Av (in km s”®) is the 
line width = 2.8 km s”®, and the geometrical parame¬ 
ter, k=126, for a density profile p oc 1/P. We there¬ 
fore find yimr= 1580 Mq for k = 126. In general, the 
clump can be stable against gravitational collapse if the 
virial ratio Msdoud/^vir 1- The virial ratio condi¬ 
tions, Msdoud/^vir < 1 and Msdoud/'^vir > 1 provide 
the signatures of unbound clump and unstable clump 
against gravity, respectively. In the present case, the 
ratio Msdoud/'^vir appears larger than unity, which in¬ 
dicates unstable clump against gravitational collapse. 
Star form ation activity is traced in this clump (see Sec¬ 


tion 3.5.2 for more details), hence the higher Psdoud value 
could be explained due to the effects of self-gravity. 

Note that the dynamical age of the H ii region is ^0.32 
Myr. Therefore, the W42 H ii region might not have in¬ 
fluenced the vicinity prior to this time period. The total 
pressure {Ptotai = ^hh + Prad + Pwind) driven by a 
massive star is found to be ^ 1.0 x 10 ”° dynes cm”°, 
which is comparable with the pressure exerted by the 
surrounding, self-gravitating molecular cloud (« 9.3 x 
10”®° dynes cm”°). This particular result provides the 
evidence that the southern condensation is not destroyed 
by the impact of the expanding H ii region. These con¬ 
clusions are valid even for the distance of 2.2 kpc to W42. 
If we compute different pressures for this distance then 
we obtain Prii = 2.1 x 10 ”° dynes cm”^, Prad = 2.3 
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X 10 dynes cm Pu,md = 1-7 x 10 dyn< 
Ptotal = 2.3 X 10“® dynes and Psdoud 

10“® dynes cm“^. 


7.8 X 


3 . 5 . Infrared excess populations 

3 . 5 . 1 . Identification and classification of infrared excess 

sources 

YSOs are identified based on their infrared excess 
emission. Here, we summarize the different schemes to 
identify and classify YSOs using MIPSGAL, IRAC, and 
WFCAM photometric data. 


IRAC-MIPSGAL bands: Guieu et al. (2010), fRe^ 


bull et al. (2011), and Samal et al. 12015) utilized a 


— [24J/[3.6] color-magnitude plot to identify YSOs using 
IRAG and MIPSGAL bands (i.e. IRAC-MIPSGAL). 
Firstly, we identified MIPSGAL 24 fim sources that 
are common with IRAC point sources. Note that since 
MIPSGAL 24 pm image is saturated near the IRAS 
location, this scheme cannot provide information of 
YSOs towards the main central ionizing area. 58 sources 
are found to be common in the IRAC-MIPSGAL bands. 
The IRAC-MIPSGAL color-magnitude diagram ([3.6] 

— [24]/[3.6]) is shown in Figure |lT] for all the identified 

sources. Using this scheme, we find 27 candidate 
YSOs (7 Class I; 3 Flat-spectrum; 17 Class II) and 31 
candidate Class III sources in our selected region. In 
Figure .0 we marked different zones occup ied by YSOs 
and C lass III sources (also see Figure 6 in Guieu et al. 
2010). Figure also shows the extinction vector with Ak 
= 5 mag which is obtained using the average extinction 
laws Ak = 0-63 2 and ki/^p,ralA k = 0.48) 


from Flaherty et al. (2007). Note that if some Class II 
YSOs sutler extinction (Ay) of about more than 20 
mag then such sources could be shifted in the location 
of the Flat-spectrum. We also checked our sources for 
possible contaminants (i.e. galaxies and disk-less stars) 
using the color- magnitude space of the SWIRE field 


(see Figure 10 in [Rebull et al.|[2011 ). 
any contam inan ts 


We do not find 
galaxies) in our selected sources 
(see Figure 11 ). One can find more details about the 


I.e. 


the work of Rebull et al 

( 12011 ). 

• I'our IRAc bands: 

dutermuth et al. 


(2009) 

tour 


devel- 

IRAC 


bands. Using IRAC colors, they also identified various 
possible contaminants (e.g. broad-line active galac¬ 
tic nuclei (AGNs), PAH-emitting galaxies, shocked 
emission blobs/knots, and PAH-emission-contaminated 
apertures). In order to identify YSOs and likely contam- 
ina nts, we adopted the vario us color criteria suggested 


by 


(2009 


Gutermuth et al. 

YS Os were further classifiei 
stages (i.e. Class I, Class 


The selected candidate 
into different evolutionary 
11, and Class III), using 


the slopes of the I RAQ SEP (ajRAn ) measured from 


3.6 to 8.0 pm (e.g., Lada et al. 
YSO classifications can also be 


20061). The details of 
ound in [Dewangan fc 


Anandarao (2011| and references therein)) The IRAC 
color-color diagram ([3.6] —[4.5] vs [5.8] —[8.0]) is shown 
in Figure [T^ for all the identified sources. Following 
this procedure, we find 39 candidate YSOs (14 Class I; 
25 Class II), 1 candidate Class III, 1285 photospheres, 
and 123 contaminants in the selected region. 


• WFCAM-IRAC bands: WFCAM-IRAC (H, K, 3.6, 
and 4.5 pm) classification method is adopted for the 
sources, that are not detected in two longer wavelengths 
of IRAC bands (5.8 and 8.0 Atm). One can find more 
details about this method in 


Gutermuth et al. 


(2009). 


In this method, the dereddened colors ([K— [3.6]]o and 
[[3.6] —[4.5]]o) w ere estimated us ing the color excess 


ratios given in [Flaherty et S) (2007). Additional 

__ t; _ F'TTn - ^ t A C _ /ill _ tt _] To cl 


conditions (i.e., [3.6]o < 14.5 mag for Class II and [3.6]o 
<15 mag for Class I) were also applied on the identified 
YSOs (Class I and Class II) to check for possible dim 
extragalactic contaminants. The dereddened 3.6 pm 
magnitudes were ob tained using observed color and the 


extinction law from Flaherty et al. (2007). We obtain 
252 candidate YSOs (16 Class I and 236 Class II) using 
WFCAM-IRAC data (see Figure [T^). 

• Three IRAC bands: One can identify additional 
protostars using only three IRAC bands (3.6, 4.5, and 
5.8 pm), when sources are not detecte d or saturated in 
8.0 pm ban d. Using three IRAC b ands, [Hartmann et al. 
(2005) and Getman et al. (2007) identified protostars 
with the criteria [3.6] —[4.5] > 0.7 and [4.5] —[5.8] > 
0.7. Following this a ppr oach, we identify 9 candidate 
protostars (see Figure 


12 :). 


H-K color excess: 'Ihose sources detected only in the 
NIR regime can be used to identify YSOs, having a large 
color excess in H—K. Red sources (having H—K > 2.35) 
were also identified using the color-magnitude (H—K/K) 
diagram (see Figure [T^). This color criterion is selected 
from the color-magnitude analysis of the nearby control 
field (size ~5l4 x 514; central coordinates: aj 2000 = 
18^38’"4U.2, 5 j 2 ooo = -06°53'53".8). We generated 
a color-magnitude (H—K/K) diagram for our selected 
control field and used all sources having detections 
in H and K bands. We obtained a color H—K value 
(i.e. ~2.35) that separates large H—K excess sources 
from the rest of the population. This cut-off condition 
leads to 185 additional deeply embedded sources in the 
complex. 

All the four schemes yield a total of 512 candidate 
YSOs in the complex. The positions of all YSOs are 
shown in Figure [T3l 

3 . 5 . 2 . Spatial distribution of YSOs 

In this section, we combine all the selected candidate 
YSOs from different schemes to examine their spatial 
distributions in the complex. To study the spatial distri¬ 
bution of candidate YSOs, we generate their surface den¬ 
sity map. The surface density map can be constructed 
by dividing the mosaic image using a regular grid and es¬ 
timating the surface density of candidate YSOs at each 
point of the grid. The formula of surface number den¬ 
sity at the z*" grid point is given by pi = {n — 1)/Ai 


(e.g. Casertano & Hut 1985[), where Ai is the surface 
area defined by the radial distance to the n = 6 nearest 
neighbor (NN). The map was created using a 5" grid at 
a distance of 3.8 kpc, which is shown as contours in Fig¬ 
ure The surface density contours are drawn at levels 
of 3cr (4 YSOs/pc^; where Icr = 1.4 YSOs/pc^), 4cr (6 
YSOs/pc^), 6 cr (8 YSOs/pc^), and 9cr (13 YSOs/pc^), in¬ 
creasing from the outer to the inner regions. More details 
o n the surface density of Y SOs can be found in the work 
of Gutermuth et al. (2009) and Dewangan & Anandarao 
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(2011). Figure 14 shows the spatial correlation between 
YSO surface density, molecular gas, and filaments. Note 
that if we use a distance of 2.2 kpc to W42, then we ob¬ 
tain the same surface density structure using the same 
contour levels (i.e. 3a, 4tT, 6cr, and 9a). However, the 
only difference is the value of la = 4.2 YSOs/pc^ at a 
distance of 2.2 kpc. 

To estimate the clustered YSO populations, we em¬ 
ployed an empirical cumulative distribution (ECD) of 
YSOs as a function of NN distance. We select a cut¬ 
off length (also referred as the distance of inflection dc) 
using the ECD, which allows to delineate the low-density 
popu l ations (see [Chava£riaetjd^ Gutermuth_et_^ 


20091 iDewangan &: Anandarao|P2Ull[ Dewangan et al. 
2015a( for more details). The analysis yielded a cut- 


off distance of dc ~0 .01262 (or ^0.862 pc) to identify 
the cluster members within the contour level of 3cr (4 
YSOs/pc^) in the entire region. A cutoff distance results 
in a clustered fraction of ~40% YSOs (i.e. 206 from a 
total of 512 YSOs). The GRS ^^CO data allowed us to 
infer the exact boundary of the W42 molecular cloud. 
In Figure [Tdl YSO clusters (i.e. gl, g2, g4, and g5) 
are spatially distributed well within the W42 molecular 
cloud, which confirms their association with W42. How¬ 
ever, the YSO cluster g3 lies close to G25.4NW. In this 
direction several molecular clouds are overlapping, and 
the molecular clump close to G25.4NW has the high¬ 
est column density of the field. Therefore, some of the 
YSOs in this direction may not be part of the cluster 
g3 and they could be associated with the clump close 
to G25.4NW. Additionally, the remaining YSO clusters 
located away from the W42 molecular cloud could be 
situated at larger distances. In general, the study of the 
distribution of molecular gas and YSOs is considered as 
a useful tool to overcome the projection effect. The po¬ 
sitions of Glass I candidate YSOs are also shown in Fig¬ 
ure Additionally, two embedded sources (positions 
and 24 p,m photometry: 1. 02000 =18:38:12.7, ^2000 = 
—6:50:01.8, m 24 = 1.79 mag; and 2. 02000 =18:38:12.3, 
^2000 = —6:52:01.5, m 24 = 4.44 mag) associated with the 
southern condensation are detected only in 24 /rm image. 
It confirms that the stars are being formed in the south¬ 
ern condensation/IRDG. The star formation activity is 
seen in all the filaments and the clumps (cNl and cSl) 
except clump cN2 (see Figures]^ and 14). The clump 
cN2 (peak N(H 2 ) ^3.5 x 10^^ cm“^) might be associated 
with the deeply embedded YSOs, which are not visible 
in the infrared regime. 

3.6. Distribution of H-band Polarization Vectors 

Figure [ 15 ^ shows the distribution of H-band polariza¬ 
tion vectors towards the molecular cloud traced by the 
integrated ^^GO intensity map. T he p olarization vectors 
of 234 stars are shown in Figure 15 1 (see Section!^ for 
more details). There are no H-band polarization detec¬ 
tions towards the dense regions (e.g. clumps cN2, and 
cSl) in the complex. The H-band polarization detections 
are observed for a few sources located within previously 
known NIR cluster including the 05-06 star in clump 
cNl. To explore the distribution of polarization data, 
we show mean polarization vectors in Figure [T^. In or¬ 
der to obtain the mean polarization data, we divided the 
spatial area into 10 x 10 equal divisions and computed 


a mean polarization value of H-band sources located in¬ 
side each specific division. The length of a vector indi¬ 
cates the degree of polarization, while the inclination of 
a vector represents the polarization equatorial position 
angle. Following the standard grain alignment mecha¬ 
nisms, the polarization vectors of background stars re¬ 
veal the s ky-projected component o f the ma gnet ic field 
direction (Davis & Greenstein 1951). Figure 16r shows 
the GPS INIR color-color diagram of sources naving H- 
band polarization detections. We find the reddened back¬ 
ground stars and/or embedded stars with (J—H) > 1.0, 
and the foreground sources with (J—H) < 1.0. Based 
on the above analysis, we find that the majority of stars 
appears behind W42 and traces the plane-of-the-sky pro¬ 
jection of the magnetic field in the W42 complex. The 
statistical distributions of the degree of polarization and 
the polarization equatorial position angles of 234 stars 
are shown in Figures [T^ and[T^, respectively. We find 
that the degree of polarization value for a large popu¬ 
lation is ^2-4%. The starlight polarimetric data show 
an ordered plane-of-the-sky component of the magnetic 
field (position angle ~15°) without additional magnetic 
field components. The overall polarization distribution is 
almost uniform towards the W42 molecular cloud. Note 
that the waist of the bipolar nebula has position angle of 
^15° (see Figure]^, which is very similar to the orienta¬ 
tion of_^e magnetic field in the plane of the sky. [Jones et| 
( 2004|) carried out K-band polarimetric observations 


towards previously known NIR cluster including massive 
star, and suggested a uniform magnetic field geometry 
threading through the entire cluster. They found a mean 
position angle of 18° and a dispersion polarization angle 
of 12° for the stars, which are located within 30" of the 
05-06 star. The distribution of GPIPS H-band polar¬ 
ization vectors towards the NIR cluster is consistent with 


the results of Jones et al. (2004). Note that the previous 
polarization study was restricted only towards the NIR 
cluster region, while the polarization study in this work 
is presented for a larger area towards the complex. 

The exact position angl e of the Galac tic magnetic field 
is not known. However, Heiles (2000) listed the opti- 
cal polarime t ry of stars iii the surroundi ng sky of W42 . 
Jones et al.j (|2004|) utilized the work of |Heiles| (|2000|), 
and suggested that the magnetic field direction in the 
diffuse ISM surrounding W42 lies in the Galactic plane 
at position angle of ^28°. 

The plane-of-the-sky projection of the magnetic field 
(i.e. mean field) associated with W42 molecular cloud is 
close to the position angle of the Galactic magnetic field. 
It indicates that the mean field is not affected by self¬ 
gravity and/or turbulence during the cloud and core for¬ 
mation processes. Gonsequently, this result suggests that 
the Galactic magnetic field appears to be the prominent 
magnetic field even within the molecular cloud (including 
embedded H ii region and filaments). Further discussion 
on the role of magnetic field is presented in Section |4^ 

On the eastern side of the bipolar lobe, we find a change 
in behavior of the polarimetric data (length and position 
angle) with respect to the polarization of stars located 
within previously known NIR cluster (see Figure 15), 
where warm dust emission is dominated. 

In order to get an idea about magnetic field strength 
in W42, we estimated the plane-of-the-sky component 


Bpos, following the equation given in Ghandrasekhar & 
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Fermi (1953): 


^pos - 1 - 




0.5 


i PI- 


( 2 ) 


-where p is the volume mass density (in g cm“^), ay is the 
gas velocity dispersion (in cm s“^), and a is the 
angular dispersion of the polarization vectors (in radi¬ 
ans). We htted a single Gaussian to the integrated 
velocity profile and obtained ay = 1.19 km s“^. In this 
calculation, we use a = 12° (mean value; as mentioned 
above) and p ~9.3 x 10“^^g cm”!^, which is compute d us- 
ing the H 2 number density (~2051 cm“^; see Section 3.41 
multiplied by molecular hydrogen’s weight (2 x 1.0(1794 
X 1.67 X 10“^^ gm) and a factor of 1.36 to account for 
helium and heavier elements. We find Bpos ~113.7 pG, 
which can be converted to magnetic pressure, Pmag (= 
i?po^(G)/87r; dynes cm“^), equal to ~5.1 x 10“^° dynes 
cm“^. If we use a distance of 2.2 kpc to W42, then we 
obtain the H 2 number density ^3541 cm“^, p ^1.6 x 
10“^°g cm“^. Epos ~148.5 pG, and Pmag "-^8.8 x 10“^° 
dynes cm“^. These values are relatively higher than that 
derived at a distance of 3.8 kpc. The values of p, Pmag, 
and H 2 number density estimated at 3.8 kpc are found 
to be about 58% of that derived at 2.2 kpc. However, 
the value of Epos calculated at 3.8 kpc is obtained to 
be about 76% of value estimated at 2.2 kpc. Note that 
we do not have polarimetric data towards the densest re¬ 
gions (i.e. cN2 and cSl clumps) in the complex, therefore 
this number can be taken as an indicative value for the 
complex. 


4. DISCUSSION 

4.1. Feedback effect of the 0-type star on the parental 
molecular cloud 


Deharveng et al. (2010) used the multi-wavelength 
data (infrared, radio, and sub-mm) to explore the trig¬ 
gered star formation at the periphery of 102 MIR bubbles 
including the bipolar nebula, N39 (W42 r egion). They 
discussed th e the oretical processes giv en in |Bodenheimer| 
et al.| (|1979|) and |Fukuda et al.| ([2000| and explained the 
presence of the bipolar morphology m N39 as due to the 
expansion of Hii region in a filament. They suggested 
that the ionized gas leaked out in two directions perpen¬ 
dicular to the dense regions. 

The analysis of CO line kinematics provided the infor¬ 
mation of the distribution of gas and its morphological 
shape in W42 complex. The position-velocity analysis 
showed the signature of the expanding H ii region (i.e. 
inverted C-like structures). The location of the massive 
star (05-06) appears to be associated with the warmest 
region (^30-36 K) compared to other surrounding re¬ 
gions (^20 and 25 K). The variation in the temperature 
inferred by the Herschel data is evident within the com¬ 
plex. A PDR region surrounding the H ii region is traced 
by the H 2 (2.12 pm) and 8 pm emissions. Two cavi¬ 
ties empty of molecular CO gas in the south-west and 
south-east directions indicated that the molecular gas is 
likely to be eroded by the ionized gas. These cavities di¬ 
rectly illustrate the interaction between the ionized gas 
of the HII region and the molecular environment in the 
complex. The W42 Hii region seems to have eroded its 
parental molecular cloud in the southern side, thus giv¬ 


ing rise to what are known as the northern and southern 
condensations. This argument is supported by the fact 
that the position angles of both the condensations being 
the same (^15°). The influence of UV photons on the 
southern condensation is also seen with the detection of 
the H 2 emission on its tip. A comparison of pressure 
contributions from different components associated with 
the massive star (i.e., Phii, Prad, and Pu,md) suggests 
that the Vhii dominates the other pressure components. 
The distrib ution of Class II c andidate YSOs (average age 
^1-3 Myr;|Evans et al.|2009|) is associated with the W42 
molecular cloud (see bection 14.3^ . The dynamical age 
for the Hii region is estimated to be ^0.32 Myr, which 
should be taken as an indicative value. This value sug¬ 
gests that a single elongated structure was present prior 
to the formation of massive star. With the time, the 
impact of the H ii region might have taken place within 
the cloud, leading to the subsequent formation of the two 
c ondensations _ 

Dale et al. (2013|) have performed the smoothed par- 


ticfe hydrodynamics (SPH) numerical simulations of the 
ionizing feedback effects from the 0-type stars on the tur¬ 
bulent star-forming clouds. They described the forma¬ 
tion of the bipolar bubble-like structures in the parsec- 
scale, due to the ionizing feedback effects. Our hndings 
also support the bipolar appearance of W42 complex as 
a result of the ionizing feedback from the 05-06 type 
star. 

In the eastern side of the bipolar bubble, a change in 
the distribution of polarization data (values and angles) 
is found. A variation in the mean angular dispersion of 
the polarization vectors (a) is also found towards the 
eastern lobe of the bipolar bubble {a ~40°) compared 
to the main molecular cloud {a ~12°). It is quite possi¬ 
ble that the ionized gas front could sweep into the dust 
grains, causing a noticeable change in the polarization 
values and angles. It seems that the ionizing gas has ex¬ 
panded in the vicinity and is responsible for the observed 
structure of the gas distribution and of the polarization 
distribution. These observational features are consistent 
with the radiation-magnetohydrodynamic simulations of 
the expansion of Hii region around an O star in a tur - 
bulent magnetized molecular cloud ( Arthur et al.||20lT ). 
The authors pointed out that the presence of magnetic 
fields is vital for the morphology of small-scale features, 
such as globules and interstellar hlaments. The results 
from the simulations also show that the expansion of H ii 
region can influence the shape of the magnetic held lines. 

All put together, our results show an imprint of the in¬ 
teraction between the H ii region and its parental molec¬ 
ular cloud. 


4.2. Role of magnetic field in W42 complex 

The knowledge of the relative orientation of the molec¬ 
ular cloud and the mean held direction allows to infer the 
role of magnetic helds in the fo rmat ion and evolution of 
the molecular cloud. In Section [3.6[ we inferred that the 
position angles of molecular cloud (northern and south¬ 
ern condensations) and the starlight polarization angles 
are consistent with the Galactic magnetic held. Note 
that the observed magnetic held (plane-of-the-sky pro¬ 
jection component) is uniform and is primarily parallel 
to the Galactic magnetic held. Therefore, it suggests the 
inhuence of the Galactic magnetic held in the evolution 
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of the molecular cloud. Similar results were also obtained 


for a massi ve star-forming region G333.6—0.2 (Fujiyoshi 
eFaI][200l| ). 


4.3. Star formation in W42 complex 

The molecular material in the region (see Figure 
could be located at different distances. However, the 
line profile along the line of sight allowed us to 
trace the molecu lar cloud associated with W42 complex. 
In Section [3.5.2[ we studied YSO clusters associated with 
W42 and also suggested the presence of a fraction of can¬ 
didate YSOs located at larger distances. The distribu¬ 
tion of clusters of YSOs (Class I and Class II) illustrates 
the star formation activity within W42 molecular cloud 
(including the filaments). In Section 3.4[ we find that the 
parental gas has been affected by tne ionized gas. The 
spatial locations of the YSO clusters in the molecular 
cloud indicate that the triggered star formation scenario 
could be applicable in W42 complex. One of the trig¬ 
gered star formation mechanisms suggests that the ex¬ 
panding H II region initiates the instability and helps in 
the collapse of a pre-existed den se clump in the molecu¬ 
lar material (e.g. Bertoldi||1989|). Gravitational instabil¬ 
ity within pre-existent condensations compressed by the 
H II region is unlikely prior to the dynamical age of the 
H II region. The Class I and Class II YSOs have an aver - 
age age of ~0.44 Myr and ~l-3 Myr (Evans et al.|2009), 
respectively. A relative comparison of these ages sug- 
gests that the age of Class I YSOs is comparable to the 
dynamical age of the Hii region (~0.32 Myr), while the 
age of Class II YSOs is higher than the dynamical age of 
the Hii region. Therefore, the evolved populations (i.e. 
Class H YSOs) are unlikely to have been the product 
of triggered formation. A small fraction of Class I can¬ 
didate YSOs is identified in the northern and southern 
condensations, that might have been influenced by the 
H II region. As mentioned before, one should consider 
the dynamical timescale of the H ii region with caution. 

We notice that the different star formation processes 
have taken place in the northern and southern conden¬ 
sations. The southern condensation harbors the densest 
(peak N(H 2 ) ^5 x 10^^ cm“^), embedded (Ay ^53.5 
mag), cold (20 K), and massive dust clump (^2564 Mq) 
in the W42 complex. The clump is associated with an 
embedded cluster of YSOs. Additionally, two embed¬ 
ded sources associated with the southern condensation 
are traced only in 24 ^m image. The NIR starlight 
polarimetric observations do not allow to obtain the 
magnetic field information within this particular dense 
clump. Therefore, we can not discuss more about the role 
of magnetic field in the star formation process. However, 
a signature of gravitational instability is ob taine d using 
the virial mass ratio analysis (see Section 3.4). Con¬ 
sequently, it seems that the cluster of YSOs associated 
with the southern condensation is most-likely originated 
by gravitational instability. 

As mentioned earlier, in the northern condensation, 
at the heart of W42 complex, there appears a parsec- 
scale cavity-like structure, which encompasses different 
early evolutionary stages of massive star formation (i.e. 
BOV star, 05-06 star, and W42-MME) (in two dimen¬ 
sional projection). The ionized cavity created by the 
UV photons of massive star(s) located in the NIR clus¬ 
ter, is not able to confine the ionized gas within it¬ 


self, which is illustrated by the observed extended ion¬ 
ized emission seen in the 20 cm map. As mentioned 
earlier, the NIR cluster contains the 05-06 star and 
two compact radio sources (i.e. G025.3824—00.1812 and 
G025.3809—00.1815). The line-of-sight velocity of the 
ionized gas in the W42 Hii region (~59.1 kms“^; Lester] 
et al.[19^ ) is very similar to the velocity of the 6.7 -CHz 
methanol maser (^58.1 kms“^; Szymczak et al. 20121 


and of the W42 molecular cloud (58-69 kms Ander¬ 


son et al. 2009). The velocity information confirms the 


physical association of molecular emi ssion, ionized emis- 


sion, and methanol maser emission. Blum et al. 
reported an extinction of Ay ~10 mag towards t 


(| 2000 | 
hiTTTR 


cluster using the average color of the brightest seven 
stars, including the spectroscopically identified 05-06 
star located within the cluster. The extinction of W42- 


MME was reported to b e Ay > 45 mag (Jones et al. 2004 
Dewangan et al. 2015b[). Based on these different values 


of extinction in W42, Jones et al. (2004) suggested the 
presence of star formation activity behind the interface 
between the H ii region and the molecular cloud. Based 
on our multi-wavelength data, W42-MME appears to be 
a more deeply embedded source compared to the O and 
B type stars. Ou r analysis is also in agreement with the 
interpretation of Jones et al. (2004). However, all these 
sources seem to be located in the same complex. 

Additionally, a “hub-filament” morphology or a fila¬ 
ment arysystem is seen in the northern condensation (see 
Section [3?^ . We find that the cavity-like structure is lo¬ 
cated at the junction of the filaments (see filaments “a, 
b, c, d, f, g, and h” in Figure]^). Note that the fil¬ 
aments have lower density compared to the hub/c lump 
associated with the cavity-like structure (see Section [3.2[ ). 
It suggests that the lower density filamentary structures 
interconnect at the high density region (see the zoomed- 
in view in Figure The spatial correlation of c lus- 
ters of YSOs and tne filaments is evident in Figure [M} 
There exist in the literature similar observational evi¬ 
dences on other cloud comple x available, such as Taurus, 
Ophiuch us, and Ros ette (e.g. Myg^[2009 Schneider et 


al. 2012|). Recently, Schneider et al. (2912) studied the 
Rosette Molecular Cloud using Herschel data and argued 
that the infrared clusters were preferentially found at the 
junction of filaments or filament mergers and their find¬ 
ings are co nsistent with the result s obtained in the sim¬ 
ulations of Dale & Bonnell (2011). The data presented 
in this work cannot throw light oii the direct application 
of this scenario in the northern condensation. A detailed 
knowledge of the motion of the molecular material along 
the filaments is required to further investigate the role of 
filaments in the formation of the YSO clusters. 

5. SUMMARY AND CONCLUSIONS 

In this paper we have studied the physical environ¬ 
ment, magnetic field, and stellar population in the W42 
complex, using multi-wavelength data obtained from the 
publicly available surveys (i.e., MAGPIS, CORNISH, 
GRS, ATLASGAL, Hi-GAL, MIPSGAL, GLIMPSE, 
UWISH2, GPS, GPIPS, ESO-VLT, and 2MASS). We 
used high resolution 5 GHz radio continuum map and 
adaptive-optics NIR images to study the small scale en¬ 
vironment of the most massive object. We utilized Her¬ 
schel temperature and column density maps as well as 
^^CO (J=l—0) line kinematics to examine the physical 
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conditions in the complex. We used different color-color 
and color-magnitude plots, as well as extinction map de¬ 
rived from the Herschel column density map, and the 
surface density analysis, to study the embedded YSOs in 
the complex. The main findings of our multi-wavelength 
analysis are the following: 

• The largest structure in W42 complex is the bipolar 
nebula with an extension of ^11 pc x 7 pc, as traced at 
wavelengths longer than 2 /im. 

• Herschel dust emissions and gas show similar spa¬ 
tial morphology with two prominent condensations (i.e. 
northern and southern) along the waist axis of the bipo¬ 
lar nebula. 

• The southern condensation is associated with the IRDC 
seen at 8.0 /xm, which is seen as a prominent bright fila¬ 
ment in emission at wavelengths longer than 70 /xm. The 
kinematics of the CO gas towards the southern conden¬ 
sation suggest that the gas is moving away at a velocity 
of ^1.5 kms“^ with respect to the center of the bipolar 
nebula. 

• The velocity of the ionized gas in the W42 H ii region 
is very similar to the velocity of the 6.7-GHz methanol 
maser and of the W42 molecular cloud, confirming their 
physical association. 

• The northern condensation at the heart of W42 com¬ 
plex contains a a parsec scale cavity-like structure which 
encompasses a BOV type object, a spectroscopically iden¬ 
tified 05-06 type object, and an infrared counterpart of 
the 6.7 GHz methanol maser (i.e. W42-MME) (in two 
dimensional projection), illustrating the presence of dif¬ 
ferent evolutionary stages of massive star formation. 


• The VLT/NACO adaptive-optics K and L images re¬ 
solved the 05-06 type star into at least three point-like 
sources within a scale of 5000 AU. 

• Two cavities of e mpty molecular gas (on scales of a few 
pc; see Figure 10) are observed in the south-west and 
south-east directions with respect to the ionizing star, 
suggesting the ionized gas has probably escaped in these 
directions. 

• The inverted C-like structures of molecular gas found 
in the position velocity maps suggest the signature of an 
expanding Hii region. 

• Herschel column density map traces two clumps in the 
northern condensation and one in the southern conden¬ 
sation. The IRDC (southern condensation) appears to 
have a peak column density of ^5 x 10^^ cm“^, which 
corresponds to a visual extinction of Ay ~53.5 mag. The 
northern condensation shows the peak column densities 
of ^3.6 X 10^^ and ^3.0 x 10^^ cm“^, which suggest 
visual extinctions of Ay ~38.5-32 mag. 

• Herschel temperature map shows a variation in tem¬ 
perature within the complex. The highest temperature 
(^36 K) is found towards the location of 5 GHz emis¬ 
sions. A PDR is traced with a temperature of ~25 K. 
The southern condensation (i.e. IRDC) appears to have 
a temperature of 20 K. 

• Parsec-scale filamentary structures are seen in the Her¬ 
schel sub-mm continuum maps which appear to be radi¬ 
ally pointed to the dense clump associated with massive 
stars, revealing a “hub-filament” system. 

• The distribution of H-band starlight polarization vec¬ 
tors shows a uniform magnetic field (plane-of-the-sky 
projection component) in the complex. The mean mag¬ 


netic field associated with W42 complex is aligned along 
the Galactic magnetic field, suggesting the influence of 
the Galactic magnetic field lines on the evolution of the 
molecular cloud. 

• 512 candidate YSOs are identified in the selected re¬ 
gion, ^40% of which are present in clusters associated 
with the molecular cloud. 

• The clusters of YSOs are distributed in the northern 
condensation including the Herschel filaments and the 
southern condensation. Additionally, the southern con¬ 
densation harbors two embedded sources that are ob¬ 
served only in 24 /xm image. 

• In the northern condensation, the cluster of YSOs in¬ 
cluding massive star is located at the junction of the fil¬ 
aments. In the southern condensation, a cluster of YSOs 
may have formed due to gravitational instability. High- 
resolution molecular line observations are necessary to 
further investigate the role of filaments in the formation 
of YSO clusters. 
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Figure 1. Herschel 70 fim emission view of W42 complex (size of the selected field ~ 15.^4 x 14.^4; central coordinates; a 2000 = 18^38^13^.9, 
<^2000 — —06°47^24'04). Contours of MAGPIS 20 cm radio emission in cyan color are superimposed with levels of 0.003, 0.006, 0.011, 0.017, 0.024, 
0.047, 0.094, 0.141, 0.189, 0.260, 0.331, 0.402, and 0.449 Jy/beam. The positions of IRAS 18355—0650 (+) and an 05-06 star (star symbo l) are 
marked in the figure. The bipolar nebula is depicted from the Spitzer 8 pm emission (dot-dashed green contour) (see [Deharveng et al.|2010[ l. The 
scale bar on the bottom shows a size of 11 pc at a distance of 3.8 kpc. The magenta box is shown as a zoomed-in view in 14gure|2^. The inset on 
the top right shows the central region in zoomed-in view, using Spitzer-lUAC 5.8 pm image in linear gray-scale (see the magenta box in the main 
figure). The 5.8 pm contour emission (dot-dashed yellow color) is shown with a level of 555 MJy/Sr. MAGPIS 20 cm emissions are also overlaid by 
cyan contours on the inset image with the same levels as shown in the main figure. The positions of a 6.7-GHz methanol maser (O) and an 05-06 
star (star symbol) are marked in the inset figure. 
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Coordinates of some referred sources in W42, as 
highlighted in Figures and[^ 
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G025.3824-00.1812 

18:38:15.3 

-06:47:52.6 

cNl condensation 
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Figure 2. a) Figure shows the central region in zoomed-in view as highlighted by a magenta box in Figure [u using GPS-K band image in 
logarithmic gray-scale. Spitzer-^.S pm contour emission (in yellow color) is shown with a level of 555 MJy/Sr. CGHNISH 5 GHz emissions (beam 
size ~iy5) are also overlaid by cyan contours on the image with the levels of 0.035 Jy/beam X (0.1, 0.2, 0.3, 0.4, 0.55, 0.7, 0.85, 0.95). The positions 
of a 6.7-GHz methanol maser emission (O) and an 05-06 star (★) are marked in the figure. The peak positions of two compact radio sources (i.e. 
G025.3809—00.1815 and G025.3824—00.1812) are shown by tri ang le symbols. A zoomed-in view of the solid magenta box is shown in Figurej^. A 
zoomed-in view of the dot-dashed green box is shown in Figure!^, b) VLT/NACO adaptive-optics and L' images (in logarithmic scale) towards 
G025.3809—00.1815. The position of an 05-06 star is highlighted by a white arrow, c) VLT/NACO adaptive-optics Kg image (in logarithmic gray 
scale) towards G025.3824—00.1812. Note that there is no VLT/NACO L image available towards this radio source. 



Figure 3. Herschel dust and molecular emissions of W42 complex, using the Hi-GAL, ATLASGAL, and GRS surveys. The panels show images 
at 160 pm, 250 pm, 350 pm, 500 pm, 870 pm, and integrated ^^CO contour map, from left to right in increasing order. In all the panels, the star 
symbol indicates the location of an 05-06 star. A zoomed-in view of the dotted box in the panel “b” is shown in Figure [^. The contour map of 
integrated ^^CO emission is shown in the velocity range of 58—69 km s~^. The CO contours are 2.37 K km s“^ X (3, 5, 7, 9, 11, 13, 14, 15). The 
magenta box in the panel “f’ is shown as a zoomed-in view in Figure [^. Based on the visual inspection of the integrated CO map, two clumps in 
the northern condensation (i.e. cNl and cN2) and one clump (cSl) in tne southern condensation are marked in panel “f’. 










Physical environment of star-forming region W42 


15 



a) -N 

46'0" 


,e-' / '' 

^ /"(g;i60) 



000) 

CO 

[c; 40)\:4 - 

CN 

^ <8' 


W ^ 


t 

50'0" 


-6°52'0" 



[250 /xm] 

1 . 1 1 1 . 1 1 1 ^Hi . 


18''38""20" 15" 10" 


a (2000) 


46'0 


47'0' 


o 

o 

o 

CM 


48' 


49' 


-6°50'0 



18'^38'^20" 15" 10" 5" 

a (2000) 


Figure 4. a) Inverted gray scale map of the waist axis of the bipolar bubble using Herschel 250 pm image. The selected area of 250 pm image is 
shown by a box in Figure!^. Filaments are shown by yellow color curves along with their designations and position angles, b) Continuum-subtracted 
H 2 image (gray-scale) at 2.12 pm of the northern condensation (size of the region ~ 4.^9 X 5.^2), as shown by a magenta box in Figure [ 3 J. The 
marked symbols are similar to those shown in Figure [^. 
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Figure 5. Spitzer 8.0 /im continuum image in log scale, as shown by a magenta box in Figure[M. The area shown in H 2 and 8.0 /xm images is 
similar. Spitzer 5.8 pm contour emission (in firebrick color) is shown with a level of 555 MJy/Sr. Filaments as seen in the Herschel maps (Figure[4^) 
are also overplotted by dot-dashed cyan color curves. The other marked symbols are similar to those shown in Figure 
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Figure 6. a) Herschel temperature map of W42 complex (see text for details). A bipolar nebula is also highlighted similar to the one shown in 
Fig urefT] A filled magen ta square refers an ATLASGAL clump, which has NH3 line parameters from the NH3 survey of dense ATLASGAL clumps 
(see |Wienen et al.|201^ . b) Herschel column density (A^(H 2 )) map of W42 complex (see text for details). The map also provides the information 
of extinction with Av — 1.07 X A(H 2 ). Two clumps in the northern condensation (cNl and cN2) and one clump (cSl) in the southern 

condensation are marked in the figure. In both the panels, the star symbol indicates the location of an 05-06 star. The Herschel filaments are also 
overplotted by solid curves in both the panels (see Figure [^). 
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Figure T. The (J=l—0) velocity channel maps of W42 complex. Each channel map is obtained by integrating five original channels between 

the velocity coverage range indicated in each panel. The color bar is shown from Icr at intervals of Icr, where the channel rms is 0.06 K km s~^. 
The position of an 0 5-06 star is marked with a filled star symbol. Two dashed curves (in black color) represent the regions empty of molecular gas 
(also see Figure[lO[l. 
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Figure 8. Two-dimensional projections of the GRS (J=l—0) data cube. Top Left: Integrated intensity map of W42 complex. The 

molecular emissions are integrated over the velocity range from 58 to 69 km s~^. The lowest gray-scale level corresponds to 20(7, with successive 
levels increasing in steps of 20a. The channel rms is 0.201 K km s~^. A zoomed-in view of the solid box is shown in bottom left panel. The 
position of an 05-06 star is marked with a star symbol. Bottom Left: Contour maps show the receding gas (65—68 kms~^; red) and approaching 
gas (59—63 kms~^; blue) with respect to that at rest (63—65 kms”^; green color) for the highlighted region in top left panel. The lowest level 
corresponds to 20a (black; la — 0.201 K kms“^), 25a (red; Icr — 0.107 K kms“^), 30(7 (green; la — 0.082 K kms“^), and 25a (blue; la — 0.123 K 
kms~^), with successive levels increasing in steps of 20<7 (black), 25a (red), 30(7 (green), and 25a (blue). The position of an 05-06 star is shown 
with a filled star symbol. Top Right: Declination-velocity map. The CO emission is integrated over the right ascension range from 18^37”^48® 
(279.45°) to 18^38^^35® (279.65°). Bottom Right: Right Ascension-velocity map. The CO emission is integrated over the declination range from 
— 06°54'00^^ ( — 6.9°) to —06°42'00^^ (—6.7°). The semi-ring-like or inverted C-like features are highlighted by dotted curves (in black color) (also 
see the text). In position-velocity maps (both right panels), the lowest gray-scale level corresponds to 10a (where la — 0.005 K degree), with 
successive levels increasing in steps of lOcr. The slices of receding, approaching, and rest gas are highlighted in red, blue, and green dot-dashed 
boxes, respectively. Cyan color dashed line shows the position of the ionizing star in the position-velocity maps. 
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Figure 9. Two-dimensional projections of the GRS (J—1—0) data cube similar to those shown in figurej^ Top Left: Integrated intensity 

map of W42 complex. The position of an 05-06 star is marked with a star symbol. Bottom Left: Contour maps show the receding gas (63.5— 
68.1 kms~^; red) and approaching gas (58.2—64.1 kms ~^; blue) with respect to that at rest (63.3—65.2 kms”^; green color). The empty molecular 
gas regions (i.e. cavities) are observed (also see figure [lO) . The position of an 05-06 star is shown with a filled star symbol. Right panels: the 
slices of receding, approaching, and rest gas are highlighted in red, blue, and green dot-dashed boxes, respectively. Cyan color dashed line shows 
the position of the ionizing star in the position-velocity maps. 
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APPENDIX 

SQL CONDITIONS TO SELECT RELIABLE POINT SOURCES FROM THE GPS CATALOG IN W42 

For all JHK: select ra, dec, jAperMagS, jAperMagSErr, hAperMagS, hAperMagSErr, k_lAperMag3, 
k_lAperMag3Err, pStar from gpsSource where ra between 279.42859 and 279.68706 and dec between -6.9096591 
AND -6.6704976 and mergedClass != 0 and (PriOrSec=0 OR PriOrSec=framesetID) AND pstar > 0.90 AND jppErrbits 

< 256 AND hppErrbits < 256 AND k_lppErrbits < 256 AND abs(jAperMag3Err) < 0.1 AND abs(hAperMag3Err) < 0.1 
AND abs(k_lAperMag3Err) < 0.1 

For only HK: select ra, dec, jAperMag3, jAperMag3Err, hAperMag3, hAperMag3Err, k_lAperMag3, 
k_lAperMag3Err, pStar FROM gpsSource where ra between 279.42859 AND 279.68706 AND dec between -6.9096591 
AND -6.6704976 and mergedClass != 0 and (PriOrSec=0 OR PriOrSec=framesetID) AND pstar > 0.90 AND hppErrbits 

< 256 AND k_lppErrbits < 256 AND abs(hAperMag3Err) < 0.1 AND abs(k_lAperMag3Err) < 0.1 
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Figure 11. Color-magnitude diagram ([3.6] — [24] vs [3.6]) of sources detected in the IRAC and MIPSGAL bands (see the text for more details). 
Class I, Class II, Flat-spectrum, and Class III sources are shown by open red circles, open blue triangles, open red diamonds, and open black squares, 
respectively. The dashed lines delineate the different spaces of these sources occupied in th e color-magnitude d iagram. The spaces of contaminated 
sources (galaxies and disk-less stars) and YSOs are separated by the dot-dashed lines_(see|R^biffletaL|201^ for more details). The arrow shows 
the extinction vector (Ajc = 5 mag) obtained using the average extinction law from [^lahert^^^ir]"|^nnY|~~' 
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Figure 12. a) Color-color diagram ([3.6] —[4.5] vs. [5.8] —[8.0]) using the Spitzer-l^KC four band detections. The “x” symbols in magenta color 
show the identified PAH-emission-contaminated apertures in the region (see the text); b) The dereddened [K — [3.6]]o vs [[3.6] — [4.5]]o color-color 
diagram using WFCAM and IRAC data; c) Color-color diagram ([3.6] —[4.5] vs. [4.5] —[5.8]) of the sources detected in three IRAC bands, except 8.0 
pm image; d) Color-magnitude diagram (H—K/K) of the sources detected in H and K bands. In all the panels, Class I and Class II candidate YSOs 
are shown by open red circles and open blue triangles, respectively. The dots in gray color represent the stars with only photospheric emissions 
(see the text for YSOs sele ction criteria). In the first three panels, the arrow shows the extinction vector obtained using the average extinction law 
from [Flaherty et al!] | |2007^ . 
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Figure 13. The positions of Class I (circles), Flat-spectrum (diamond), and Class II ( triangles) candidate YSOs identified within our selected 
region. The candidate YSOs identified usin g th e IRAC-MIPSGAL scheme (see Figure [ll^ are shown by magenta color, while the candidate YSOs 
selected using the other schemes (see Figure [l2[ l are highlighted in blue (Class II) and red (Class I) colors. The background map shows an inverted 
gray scale Herschel 70 /^m continuum image. The position of IRAS 18355—0650 source is shown by a “-h” symbol. A bipolar nebula is also 
highlighted similar to the one shown in Figure 
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Figure 14. Overlay of YSOs surface density contours on the integrated (J=l—0) emission contour map. The surface density contours are 

drawn at 3(T (4 YSOs/pc^, where Icr—1.4 YSOs/pc^), 4cr (6 YSOs/pc^), htr (8 YSOs/pc^), and Qct (13 YSOs/pc^), increasing from the outer to 
the inner regions. The background map is similar to the one shown in Figure[^. Class I candidate YSOs are shown by open red circles. Five YSO 
clusters are labeled in the main molecular cloud. Two sources are marked by white asterisk symbols, which are detected only in 24 fim map. The 
Herschel filaments are overplotted by orange color curves (see Figure The position of an 05-06 star is marked with a filled star symbol. 
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Figure 15. a) Overlay of GPIPS H-band polarization vectors (in red color) of 234 stars with UF = 1 and P/c7p > 2 on the integrated 
emission map. The background map is similar to the one shown in Figure Vectors are drawn for the region marked by a dot-dashed box in 
figure, which are retrieved from the GPIPS fields GP0612 (/ — 25°.319, b — —0°.240) and GP0626 (/ — 25°.447, b = —0.°165). The length of each 
vector represents the degree of polarization (a reference vector of 5% is shown in the lower left corner of the map). The orientations of the vectors 
represent the equatorial position angles of polarization, which show that the surrounding interstellar magnetic field lies in the Galactic plane at 
position angle of 28°. b) Overlay of mean polarization vectors on the inverted Herschel 250 gray scale map. The polarization spatial area is shown 
by a box in Figure [T^. The mean polarization data are obtained by dividing the polarization spatial area into 10 X 10 equal divisions and, a 
mean polarization value of H-band sources is estimated inside each specific division. A bipolar nebula is also highlighted similar to the one shown 
in Figure In both the panels, the star symbol indicates the location of an 05-06 star. 
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Figure 16. a) GPS NIR color-color diagram (H—K vs J—H) of the sources hav ing H-band polarizatio n observations. The solid curves represent 
the unreddened locus of main sequence stars (cyan) and giants (magenta) (from|Besse^^^Brett|1988jl. The long-dashed straight lines show the 
extinction vectors with Ak — 3 mag. The extinction vectors are derived from [Colier^^^iri't^dlll extinction laws {Aj/A k — 2.94 and Ah/Ak 
— 1.72 for California Institute of Technology (CIT) system). The color-color diagram is divided into three differe nt subregions, nam ely “I”, “II”, 
and “III”. The YSOs are identified in “I” and “11” subregions. Classical T Tauri (CTTS) locus (in CIT system) | |Meyer et al.|l997| is plotted as 
a solid yellow line. The loci of unreddened dwarf (Bessell fc Brett (BB) system), giant (BB-system), and colors are converted into CIT system 
using transformation equations provided by ICarpenterl l|2001|l. b) The histograms of degree of polarization; and c) equatorial position angles of 
polarization obtained for 234 sources as mentioned in Figure 1^. 




























